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ABSTRACT 

New  phenomena  are  observed  theoretically  in  nonlinear  optical  responses  of 
polvdiacetylene- toluene  sulfonate  (PTS)  irradiated  by  a  strong  laser  beam.  Within  a 
two-level  model  of  PTS,  the  transient  behavior  of  the  induced  susceptibility  is 
investigated.  Optical  nutation  is  found,  and  the  optical  Stark  blue  shift  and 
bleaching  are  also  found  in  qualitative  agreement  with  experiments.  In  the  steady 
state,  a  new  type  of  optical  tristability  mediated  by  phonons  is  found.  A  possible 
mechanism  responsible  for  this  tristability  is  discussed. 

1 .  INTRODUCTION 

Nonlinear  optical  processes  have  been  employed  to  realize  light  frequency 
conversion,  light  modulation  and  switching,  optical  logic  and  memory  storage  in 
optical  computing,  signal  processing  and  image  analysis.1  Because  of  such  immense 
potential  applications,  it  has  been  one  of  the  most  active  areas  of  research  for 
both  basic  understanding  of  the  mechanism  and  technological  development  of  electro- 
optical  or  all-optical  devices  for  communication  and  computer  systems. 

Organic  polymers  are  known  to  exhibit  giant  nonlinear  optical  susceptibilities 
with  fast  response  time.2  In  particular,  the  third-order  susceptibility  *  '  of 
conjugated  polymers  has  attracted  much  attention.  It  is  expected  to  plav  an 
important  role  in  phase  conjugation,  optical  bistability  and  all-optical  modulation. 

The  large  optical  nonlinearities  in  conjugated  polymers,  first  predicted  over 
ten  years  ago,3  have  been  observed  in  numerous  systems  in  recent  years.  The  largest 
X'  is  recorded  for  trans -polyacetylene4  (CH)  for  which  a  number  of  theoretical 
investigations  have  h.een  carried  out  very  recently.  On  the  basis  of  a  one - 
dimensional  mcdel,5,6  x  has  been  calculated  for  (CH)  by  means  of  the  Keldysh 
Green's  function,7  and  fo^  pristine  (CH)^  with  neutral  s'olito^s.8  Furthermore,  the 
chain-length  effects  on  x  for  (CH)  have  also  been  studied. 

Meanwhile,  nonlinear  optical  responses  in  polydiacetylene  (PDA)  have  attracted 
much  attention.  The  nonlinear  refraction  index  n.  of  PTS  films10  and  PDA-4BCMU11 
have  been  measured,  and  phonon- mediated  excitonic  processes  in  PTS  have  been 
observed  in  transient  nonlinear  optical  response12,13  and  in  an  optical  Stark  effect 
experiment . 1 4 

We  report,  in  this  paper,  theoretical  studies  of  the  transient  nonlinear 
optical  response  of  virtual  excitons  in  PTS  and  a  new  type  of  optical  tristabilitv 
of  PTS  in  an  optical  cavity. 

H 


□  □  C_L 


2.  TRANSIENT  OPTICAL  RESPONSE  OF  PTS 


The  investigation  of  coherent  nonlinear  optical  processes  caused  by 
semiconductor  excitations  in  the  transparent  region  well  below  the  absorption 
edge15-18  has  been  extended  to  polymers.12-14  It  is  known  that  virtual  excitons  are 
created  when  the  frequency  of  an  incident  laser  beam  is  tuned  below  the  exciton 
resonance.  Such  excitonic  effects  can  be  explored  more  naturally  by  probing  the 
susceptibility18  md  looking  for  transient  changes  that  persist  as  long  as  the 
excitation.  Proviued  that  the  sample  is  strongly  excited  well  below  the  conduction 
band  edge,  the  real  part  of  the  susceptibility  can  be  examined  by  measuring  the 
index  of  refraction,  and  the  imaginary  part  is  explored  by  the  absorption  spectrum. 

(3) 

As  PDA  (and,  in  turn,  PTS)  possesses  fairly  large  x  and  small  transmission 
loss  a,  which  can  be  reduced  to  as  low  as19  a  <  1  db/cm,  PDA  can  have  a  rather  large 
ratio  x  -  In  addition,  PDA  is  relatively  easy  to  use  in  the  construction  of 
waveguides.  In  this  section  we  take  PTS  as  an  example  and  study  the  transient 
behavior  of  its  induced  optical  susceptibility  due  to  irradiation  by  a  strong 
pumping  laser  and  a  weak  probing  laser.  The  two  lasers  may  have  pulse  duration  of 
femtoseconds  to  picoseconds.  For  the  PTS  system  we  adopt  a  two-level  model  to 
describe  its  ground  state  and  exciton  state.  The  lifetime  of  the  exciton  is  also 
introduced  phenomenologically. 

2.1  Two- level  model 


The  two  states  of  PTS  are  denoted  by  | +>  and 
transition  and  define  the  dipole  operators 

S'  -  |-X+|  , 

S+  -  l+x-l  . 

The  inversion  operator  is  given  by 


|  (l+x+l 


-X- 
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->.  We  assume  a  dipole 

(la) 

(lb) 

(2) 
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For  simplicity,  the  effect  of  phonon  mediation  “  is  ignored  here.  (Its  inclusion 
will  be  discussed  and  published  elsewhere.)  Thus  we  consider  a  single  virtual 
exciton  and  write  the  Hamiltonian  in  the  rotating-wave  approximation  as 


H  -  H  +  H' 
o 

,  iu  t 

Ho  -  J  <„oEpe  V  S-  +  h.c.) 

1  1Wtt  - 

H'  “  -  2  (M0Ete  S  +  h.c.)  , 


(3a) 

(3b) 

(3c) 


where  is  the  exciton  frequency,  n  is  the  u  vole  matrix  element,  and  E  and 
cj  areX  the  amplitude  and  frequency  of  the  pump  and  test  fields,  respectively. 
Hite  we  have  neglected  the  momentum  dependence  of  virtual  excitons.13,23  This 
Hamiltonian  uses  dipole  or  projection  operators  to  describe  the  dipole  transition 
between  the  ground  state  and  the  excitonic  state  of  the  material  by  external  fields. 


It  is  different  from  that  of  Ref.  13,  where  creation  and  annihilation  operators  were 
used  for  virtual  excitons. 


We  assume  further  that  initially  the  system  is  in  its  ground  state  with  no 
virtual  exciton  present.  Thus  we  have  |^(0)>  -  |->  at  t  -  0.  For  any  later  time  t 
>  0,  the  state  can  be  represented  by 

|0(t)>  -  C+(t) |+>  +  C_(t) | ->  .  (4) 

Plugging  (4)  into  the  Schrodinger  equation  with  the  Hamiltonian  (3),  we  find  the 
equations  satisfied  by  the  amplitudes  C+  in  the  rotating  frame  as 

■ii  t  - iA  t 

C+(t)  -  in*  e  P  C  +  in*  e  C  C  -  7C+  (5a) 

iA  t  iA  t 

C  (t)  -  in  e  P  C+  +  int  e  t  C+  ,  (5b) 

where  we  have  defined  n  —  4QE  and  A  —  u>  -  u>  .  In  addition,  a  damping 

rate  7  has  been  introBhced  toPdescribe  tRfeCver  ^’sh  ortX lifetime  of  the  virtual 
excitations . 


Since  the  test  field  is  in  most  cases  much  weaker  than  the  pump  field,  we  can 
safely  apply  perturbation  theory  to  solve  (5).  The  zeroth-order  equations  are 
obtained  by  setting  E  =  0  in  (5),  namely, 
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In  what  follows,  we  assume  for  simplicity  p  —  p*  and  n 
conditions  C°(0)  -  0  and  C°(0)  =  1,  we  find°the  solutioR^ 
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(7a) 

(7b) 


where  we  have  defined  the  complex  Rabi  frequency 

0  -  2  KAp  +  iy)2  +  4n2]H  .  (8) 

The  real  part  of  the  Rabi  frequency  is  given  by 

nR  -  [ (7x2  +  y2  +  x)/8 ] ^  ,  (9a) 

which  determines  the  oscillation  frequency,  and  the  imaginary  part  by 

n.  -  - [ (7x2  +  y2  -  x)/8]h  ,  (9b) 

2 

which  determines  the  envelope  of  the  oscillation.  In  (9)  we  have  defined  x  -  A  + 

4n  -  7  and  y  -  2A  7,  with  A  <  0.  p 

P  P  P 


The  solutions  up  to  first-order  perturbation  can  be  obtained  by  inserting  the 
zeroth-order  solutions  (7)  into  the  right-hand  side  of  (5).  The  results  are 


C+(t)  -  C+(t)  +  SC+(t)  ±  a4 


(10) 
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The  integration  constants  A,  in  (10)  are  determined  by  the  initial  conditions. 


The  expectation  value  of  the  dipole  moment  of  a  virtual  exciton  is 

-  iw  t  iw  t 

fi(  t)-u<S'(t)>e  X  -  M  C*C  e  X 
o  o  -  + 

*  *  - iw  t 

-  ^  (C°  +  SC  +  A*)(C°  +  SC  +  At)  e  X  (12) 

O  -  -4*  +4* 

It  is  evident  that  the  dipole  moment  induced  by  the  test  beam  is  given  by 


-  iw  t 

Sp(t)  -  a»  (C°*5C  +  C°5C*  ) e  X 
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(13) 


Therefore,  the  nonlinear  optical  susceptibility  experienced  by  the  test  beam  i 
simply 
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Z1  “  ^iAp  '  '  iAtJ2  +  r<2 

Z2  “  [2(iAp  *  7)  ’  iAt]2  +  n*2  •  (15b) 

We  have  introduced  in  (14)  the  optically- induced  transient  virtual  exciton  density  n 
which  may  be  assumed  to  obey  the  rate  equation23 

"  =•  r/a  Ip/(Mw  )  -  yn  ,  (16) 

where  p  denotes  the  quantum  efficiency  of  virtual  exciton  generation  per  absorbed 
photon,  I  -  E  c  n  /8n  is  the  pumping  intensity  with  n  as  the  linear  refractive 
index,  anR  a  rlpresents  the  absorption  coefficient  whict?  is  assumed  to  be  a  constant 
in  the  present  case.  In  practice,  however,  a  may  be  related  to  the  density  n  to 
some  extent. 


Equation  (16)  can  be  solved  by  a  solution  of  the  form 
n  ”  ns<'1  "  e  1  (17) 


wnere  n  —  pal  /).  Once  again,  the  initial  condition 

assumed”.  The^  tratPs ient  behavior  of  the  nonlinear  optical 
calculated  numerically  as  a  function  of  time  and  the  detuning  A 
presented  and  discussed  in  the  next  section. 


n(0)  =  0  is  already 
susceptibility  \  Is 
.  The  results  are 


2.2  Optical  nutation  and  hole  burning 


In  our  numerical  study  we  consider  PTS  and  use  the  unit  jrf  -  1.  Other 
^  a  r  aji  eters^are  chosen  as12-14  A  =  0.1  eV ,  0  =  0.02  eV  ,  y  =  0.01  eV, 
2 1  A*  |  eV  and  »  2  eV.  ifte  real  and  imaginary  parts  of  x  are  computed 
separately.  figure  I  shows  snapshots  at  different  instants  of  time  ofCthe  variation 
of  x  t  with  the  detuning  A^_  in  the  unit  of  eV.  The  variation  with  time  for  fixed 
detuning  is  shown  in  Fig.  2. 

It  is  clearly  seen  from  Fig.  1  that  both  the  differential  absorption  and 
refraction  are  direct  effects  of  the  virtual  excitons  induced  by  the  pumping  field. 
Nonzero  susceptibility  starts  to  show  up  as  |A  j  -*  0.1  and  approaches  its  maximum  as 
A  -»  0 ,  where  the  test  field  is  at  resonance  with  the  exciton.  This  indicates  that 
t.Rere  is  energy  transfer  from  the  pump  field  to  the  test  field  via  virtual  excitons, 
in  qualitative  agreement  with  experiments.12,13 

Optical  nutation  can  be  observed  in  polymers  in  the  ultrashort  time  regime 
when  virtual  excitons  are  generated  by  a  pump  field  tuned  below  the  exciton 
resonance.  This  situation  can  be  seen  from  Fig.  2.  It  is  also  clearly  shown  that 
the  oscillation  centers  of  both  Rext  and  Imxt  change  with  time.  This  is  quite 
different  from  nutation  phenomena  predicted  in  Refs.  21  and  23  for  semiconductors, 
where  the  oscillation  amplitude  decreases  monotonically  like  damped  oscillators. 

Since  it  takes  time  for  the  pump  field  to  build  up  enough  virtual  excitons  and 
since  the  virtual  excitons  are  short-lived,  their  effects  are  appreciable  and  stable 
most  of  the  time  during  the  pulse  (on  the  order  of  picoseconds).  The  response  time 
can  be  anywhere  between  10  —  100  femtoseconds  according  to  our  calculation.  Both 


the  cbsorotion  and  refraction  parts  show  up  almost  right  awav  and  taper  off 
rather  slowly  after  reaching  their  peaks. 

•.ddition,  we  can  also  observe  from  these  results  that  increasing  the 
itens:  if  the  pump  field  leads  to  a  deeper  and  wider  hole  in  the  absorption 

spectr .  which  is  qualitatively  in  agreement  with  the  steady-state  results  of 

Ref.  10.  What  seems  to  be  more  interesting  is  that  we  find  a  blue  shift  of 
0.005  eV  in  our  calculation,  even  though  we  have  not  considered  the  effect  of 
phonons  .  A  more  careful  treatment  with  phonon  mediation  included  can  be  found 
in  Ref.  22.  This  is  the  electronic  state  optical  Stark  shift  observed  in  a 
recent  experiment14  in  which  the  differential  transmission  of  polydiacetvlene- 
toluene  sulfonate  has  been  measured  by  means  of  coherent  inverse  Raman 
spectroscopy.  It  is  demonstrated  in  Ref.  14  that  this  ac  Stark  effect  is 
important  in  the  determination  of  the  nonlinear  optical  response  of  PTS  even  in 
che  small  signal  limit.  The  effect  is  reproduced  theoretically  in  Ref.  14  by 
including  an  extra  term  to  modify  the  exciton  resonance  frequency  in  the  simple 
model  of  Schmitt-Rink  et  al.13  However,  the  treatment  in  both  Refs.  13  and  14 

assumes  a  steady-state  solution,  while  the  experimental  conditions  are  more 

transient- like . 

3.  PHONON-MEDIATED  OPTICAL  TRISTABILITY  IN  PTS 

We  now  turn  our  attention  to  a  special  case  in  which  a  sample  of  PTS  is 
placed  inside  an  optical  cavity,  and  investigate  its  response  to  a  laser  beam 
directed  into  the  cavity.  In  particular,  we  look  for  the  dependence  of  the 
output  field  intensity  on  the  intensity  of  the  incident  laser  field.  To 
describe  our  system,  we  assume  the  model  Hamiltonian  of  Ref.  13  with  an  extra 
term  representing  the  cavity  photon-exciton  interaction.  Thus 

H  -  Qa^a  +  b^b  +  u>  c^c  +  Aclc(b  +  b^) 
ph  ex 

+  ig(a^c  -  c^a)  +  i(a^£e  -  a£*e^C)  ,  (18) 

where  a^(a),  b^(b)  and  c^(c)  are  the  creation  (annihilation)  operators  for  the 

cavity  field,  phonon  and  exciton,  respectivly,  with  corresponding  frequencies  0 , 
•j>  ^  and  w  .  The  phonon  -  e  xc  i  ton  coupling  constant  is  denoted  by  X,  and  the 
exciton- f ie¥d  coupling  constant  by  g.  The  coherent  driving  field  is  assumed  to 
have  the  same  frequency  Q,  and  its  amplitude  is  E  whose  spatial  dependence  is 
ignored  completely.  For  simplicity,  we  have  included  only  one  (  f  the  phonon 
modes24  which  couples  mos*-  strongly  to  the  exciton,  with  momentum  dependence 
neglected. 1 3 

In  the  usual  zero  -  fluctuation  limit,  quantum  mechanical  operators  can  be 
replaced  by  their  expectation  values  a  —  <a>,  fi  —  <b>,  rj  -  <c>  and  n  —  <c  c> , 
and  in  the  limit  of  high  Q  cavity,  the  phonon  and  exciton  variables  are 
eliminated  adiabatically . 2 5  Without  going  into  the  details  which  will  be 
published  elsewhere,28  we  discuss  only  briefly  the  results  here.  The  equation 
of  motion  for  the  field  variable,  expressed  in  the  unit  such  that  7  -  1.  takes 

the  form 

-  '2g2  E  ■  «cf(|»c!2)  .  (19) 

where  the  function  f  is  defined  as 


f(x)  -  1  + 


(20) 


1  -  i [A  -  X  n^(x) ] 


1  +  f  A 


A  n  ( x )  j ' 
P  c 


In  Eq .  (19),  we  have  introduced  for  convenience  the  phenomenological  damping 

rates  r.  ,  y  and  7  and  dimensionless  variables  r  —  k  t  ,  E  =  E/k  and  a  - 


d  n  ex 

Thfe  field  intensity  is  expressed  as  I  = 


and  the  detuning  as 


In  the  steady  state,  we  set  a  —  a  -  constant  in  (19)  and  (20).  The 
^  c  s 

?adv-state  cavity  field  intensity  can  then  be  found  straightforwardly  as 


1  +  Y  (I.  ) 

I  = - - - — -  I 

s  2.2  2  in 

(1  +  g  A)  f  Ys(i.n) 

where  we  have  defined  the0driving  field  intensity  1^ 


field  intensity  I.  -  2 g^ 1 E I ^  and  Y  (I.  )  -  A 

J  in  oil  sin 

It  can  be  shown  from  (21)  that  the  steaav-state 


-  X  n  ( I  .  )  w  i  th  X 

.ps.  in . 
virtual  exciton  n 


2X  /u 


latisf ies  Pthe  equation 


2 ( A/ X  ) n  + 
P  s 


2  2  2  -’ 

I  +  (1  +  g  /*)  ; X  n 

P  S 


I.  X 
in  p 


’.•.'hen  the  conditions  are  such  that  Eq .  (22)  has  three  distinct  roots,  there  will 

he  three  output  I  for  each  given  I.  .  A<-  has  been  discussed  in  Ref.  26,  the 

r  s  0  in 

conditions  of  multi-real  roots  of  (22)  require  that  tristabilitv  occurs  only 
when  the  cavity  field  is  tuned  sufficiently  below  the  exciton  resonance  and  when 
-here  exists  phonon-exciton  coupling.  This  clearly  shows  that  the  tristabilty 
is  mediated  by  phonons. 


Numerical  calculation  of  Eq .  (21)  has  been  performed  for  PTS  for 

which13,14  u  -  5.16  ari^l  X  -  2.0.  The  results  are  plotted  in  Fig.  3  for 

various  choices  of  A  and  g V* .  It  is  found  that  all  the  multiple  solutions  are 

stable  except  for  those  around  the  turning  points  A  and  B  in  Figs.  3(b)  and 

3 1  c )  .  As  I.  increases  continuously,  passes  the  turning  point  and  then 

decreases  ,  £Re  output  I  will  follow  a  reversed  hysteresis  loop.  The  middle 

branch  can  onlv  be  reached  if  I.  reverses  its  course  of  change  at  the  turning 

in  0 

roint. 


It  is  emphasized  that  the  dispersion  effect  plays  a  crucial  role  in  this 

kind  of  tristabilitv.  When  I.  increses  bevond  a  certain  level,  n  becomes 

-  in  ■  s 

large  enough  that  the  lattice  relaxation  is  sufficiently  high  for  the  refractive 

index  to  deviate  from  the  cavity  resonance.  Consequently,  the  cavity  field 

switches  to  the  lower  branch.  On  the  other  hand,  such  tristabilitv  cannot  occur 

when  the  Q-value  of  the  optical  cavity  is  excessively  high.  As  a  final  remark, 

we  remark  that  this  tristability  is  found  by  considering  only  one  phonon  mode  in 

PTS.  It  is  therefore  expected  to  find  optical  multistability  if  more  phonon 

modes  are  included. 
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Figure  3.  Intensity  of  the  cavity  field  I  versus  the  input  field 

2  S  2 

intensity  I.  :  (a)  A  -  6,  g  /k  -  0.1;  (b)  A  6 ,  g  /«  -  0.5;  (c)  A-4, 

o  1-0  2 

-  0.5;  (d)  A  -  5,  g  /k  -  2 .  The  dashed  lines  in  (b)  and  (c)  are 

hypothetical  and  drawn  onlv  to  complete  the  reversed  hysteresis  loop. 
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